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Abstract 

We remark that the weak coupling regime of the stochastic sta- 
bilization of 2D quantum gravity has a unique perturbative vacuum, 
which does not support instanton configurations. By means of Monte 
Carlo simulations we show that the nonperturbative vacuum is also 
confined in one potential well. Nonperturbative effects can be assessed 
in the loop equation. This can be derived from the Ward identities of 
the stabilized model and is shown to be modified by nonperturbative 
terms. 
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1. Introduction 

It is well know that discretized 2-D quantum gravity coupled to conformal 
matter with c < 1 can be represented by a zero-dimensional matrix model 
[H, H|, [3| and that the expansion in powers of 1/N of the matrix model, where 
N is the dimension of the hermitian matrix field, defines the topological 
expansion of the original discretized 2-D quantum gravity The partition 
function of zero-dimensional matrix model has the form 

Z = J d$ exp{-NTrW ($)} (1) 

where W is the potential and $ is an hermitian matrix field. The continuum 
limit of the discretized 2-D quantum gravity is achieved by taking the limit 
of the coupling constants of the potential W to some critical set of coupling 
constants for which the partition function (Jl|) has non-analytic behaviour. 
For pure gravity the argument of the integral (|1|) is an increasing function for 
large values of $ and the integral does not exist for finite N, but in the large 
N limit the integral exists up to some value of the coupling constant. This 
is the critical value of the coupling constant, and the non-analytic behaviour 
developed makes the partition function above it not well-defined. The matrix 
model is defined only in the large N limit and its expansion in powers oil/N 
represents the topological expansion of the quantum gravity, but the matrix 
model does not define uniquely the non perturbative 2-D quantum gravity 

I- 

Stochastic stabilization || § provides, on the other hand, a way of map- 
ping the model (|l|) into one which is defined for all values of iV and the 
coupling constant, while reproducing the perturbative expansion in powers of 
1/N. There has been much hope that this stabilization of the original matrix 
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model could lead to an unambiguous nonperturbative definition of quantum 
gravity. If these expectations have not been met, it has been in part because 
of the poor understanding of nonperturbative effects in the stabilized the- 
ory. As we review afterwards, this is defined in one more dimension than the 
original matrix model, setting therefore the problem into the framework of 
quantum mechanics. In fact, in the zero-dimensional model ([!]) the source of 
nonperturbative effects lies on the transmission of eigenvalues through the 
walls of the confining potential well J7j, so that one could expect the same 
kind of phenomenon in the stabilized theory. The main point of this article, 
however, is that the ground state of the one-dimensional model does not bear 
tunnelling between different wells of the potential. We are able to prove this, 
for the simplest matrix model, in perturbation theory as well as in the non- 
perturbative regime (using a Monte Carlo approach in this latter case). This 
leaves open the proper interpretation of the nonperturbative effects. Their 
significance can be assessed otherwise, since, as also shown in the paper, they 
modify the form of the loop equation in the stabilized theory. 
2. Perturbative approach 

The stochastic stabilization of the matrix model introduces a positive 
definite hamiltonian 



This hamiltonian is well defined for all values of N and the coupling constants. 
The zero mode of the stabilized theory is given by 
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its norm being the partition function of the original matrix model. Hence, 
the original matrix model is well defined if and only if this zero energy state 
is a normalizable state. When this is the case, corresponding obsevables in 
both theories coincide 

(Q)stab = jQ\V\ 2 d$ = ^J Qexp{-iW($)}d$ = (Q) matrix (4) 

We can use the stabilized theory to calculate the observables of the origi- 
nal matrix model, with the ground state energy playing the role of an order 
parameter which gives us information about the range of definition of the in- 
tegral ([!]) • The non-analytic behaviour of the original matrix model becomes 
now non-analytic behaviour in the ground state energy of the hamiltonian 
(HD- For pure gravity in the large N limit, below the critical point the zero 
energy state is a normalizable state and the ground state energy is zero, 
while above the critical point the zero-dimensional matrix model does not 
exist and the ground state energy of the stabilized hamiltonian is greater 
than zero. The stabilized theory exists anyway for all values of N and the 
coupling constant, even if the original matrix model is ill defined. 

We spend a few time showing that, despite the strange critical expo- 
nent for the leading contribution ||, the topological expansion of the ground 
state energy can be organized in terms of the scaling variable of the zero- 
dimensional matrix model. Dealing with the 1/N expansion we will also show 
our main conclusion in the weak coupling regime, i.e. that the eigenvalues 
of $ in the ground state are all confined in the same well of the potential in 
(0). Let us consider the cubic matrix model given by 

W($) = Tr$ 2 - — Tr$ 3 , (5) 
3 
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where $ is the N dimensional hermitian matrix. We consider throughout this 
section the large N limit of the model, for which (ffl) is still defined within 
a certain range at g > 0. It is well-known that the one-dimensional matrix 
model is mapped into a gas of free fermions ||, which are the eigenvalues {A n } 
of the $ matrix variable]]. The Fokker-Planck hamiltonian of the stabilized 
model (H) becomes the sum of N one-particle hamiltonians || ||, [12], [13], |~| 

K = ~\jp^2 + \{9 2 K, - 2gX 3 n + % + 2^A n - 1} (6) 

The one-particle hamiltonian (Mj has a main well and a local minimun below 



the critical coupling constant g c = y l/(6\/3), and only one well above it (see 
figure 1). 

The ground state energy of the Fokker-Planck hamiltonian is 

N-l 

E=J2 e n (7) 
n=0 

where {e n }o< ra <Ar_i are the first iV eigenvalues of the one particle hamiltonian 
(0). In the large N limit we can write down for it an integral representation, 



taking into account 1/N effects. From the Euler-Maclaurin sum-formula [15 
the total energy is 

JV_1 r E F 1 1 11 11 

E=Ye n = dEEp(E)--E* F + -E*n + ; ; — - + ••• (8) 

n =o J** Py ' 2 F 2 12 p(E* F ) 12 p(E* ) 1 ; 

where Eq is the first eigenvalue eo, E* F is the value of the analytic continuation 
of e n from n < N — 1 to n = N, and p(E) is the density of states. 



1 In other matrix models, the hamiltonian of the stabilized theory has interaction terms, 
but one can perform a Hartree-Fock approach to obtain the exact ground state energy in 
terms of suitable one-particle states pi ITol, O . 
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From (^j) we see that 1/N play the role of h, and the large N limit is a 
semiclassical expansion. The density of states is given by 



P(E) = ^N(E) 
The distribution function of states is given by [[l 

N r b — — — i 



(9) 



b „d 2 V 
dX 



dX 2 



2(E - V) 



1 

N 3 



(10) 



a and 6 being the classical turning points. Ep and Eq obey the quantization 
conditions 

HK) = \ (11) 

and 

N(E F )=N+l (12) 



where N(E) is the distribution function N(E) restricted to the main well 
of the stabilized potential. In general, one has to keep track of higher order 
terms in the sum formula @ altogether with the perturbative expansion 
(|10D, in order to obtain correctly the 1/N 2 expansion of the total energy. 
From (H) and (10) the total energy is 



NE = N 2 ( E + E 2 -- + 0[ — 



N 2 



N 4 



where 



£; = 4 0) -^-/ 6 dA[2(4 0) -y)] 

D7T J a 



Eo 



24tt 



b „d 2 V 
dX 



rfA 2 



2(4 0) 



V) 



1 

24 



1 r b 
- / dX 



7T Ja 



2(4 0) 



V) 



(13) 

(14) 
(15) 
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Here Ep is the Fermi energy in the planar limit and is given by the condition 

- [ b d\^2(EP-V) = l (16) 

77 J a 

We expect the total energy to be zero for g < g c to all perturbative 
orders f|, which we have actually checked to second order ml/N expansion. 
The leading order of the total energy ( |14) is zero up to the critical coupling 



g c = yl/(6\/3), and greater than zero above it. The critical behaviour of 
this quantity agrees with the exponent found in the cuartic model || 

E ~(g-g c ) n / A g>g c (17) 

Regarding the subleading contribution fll5|) we find the critical behaviour 

E 2 = g<g c (18) 
E 2 ~ (9-9c) 1/4 9>9c (19) 

The appropriate double scaling limit g — > g c and N — > oo is such that 
z = N(g c — g)i remains finite 0. We see, in fact, that the topological 
expansion of the total energy organizes above the critical point as a power 
series in the scaling variable z of the original matrix model 

NE = N 2 (g - g c )^ (b 1 + B 2 1 +■■) (20) 
V N\g-g c )2 J 

The striking result, however, concerns the position of the Fermi level 

Ep. To leading order of the 1/N expansion, Ep is characterized by the 

condition (|16"D , which places it precisely at the level of the local minimum of 

the potential, all the way up to g c JIT). The first order correction E^ can 

be computed from the quantization condition 

^/>yi(^ +0 (i)= n+ i (2D 



bearing in mind that we have to fill the first iV levels from eo to e^-i as in 
(0). We find then 

1 



E 



(i) 

F 



2 



l r b l 
- dX- 



(22) 



n Ja p(Ef - V) 
which is a negative quantity. This means that, at least in the weak coupling 
regime, the ground state appears to be confined to the central well of the 
potential, and that there is no issue about tunnelling to the region around 
the local minimum. One may think that this only points at the unfeasibility 
of discussing nonperturbative effects in the very framework of the 1/N ex- 
pansion. Actually, this approximation implies taking the limit N — > oo from 
the start, so that it would be conceivable that a more sophisticated way of 
tuning the double scaling could lead to a different physical picture. In the 
last section we will return to the question of the localization of the ground 
state of the model at finite N. 

3. Observables and the loop equation 

Below the critical point, the observables of the stabilized theory and the 
matrix model have to be the same, and may be calculated as follows Jl2|, [11] . 

We add to the potential V an auxiliary term f5X n and, using the Hellmann- 
Feynman theorem, write the following formula for observables of the stabi- 
lized theory 

= . (23) 



N ' \ dp J 



j3=0 



Taking into acount (|I3"D, (|23| ) and results from appendix A, we perform 
the calculation of the observables of the theory up to second order in j± . To 
first order all observables are given by 

(tr$ n ) = N 2 f b d\X n p(\) (24) 

J a 
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and 

p(X) = -^2{E%-V) (25) 

is the fermionic density and is equal to the eigenvalue density of the original 
matrix model in the planar limit for g < g c fpTOf] . Now we perform the 
limit g — > g c , which represents the continuum limit of the discretized 2-D 
quantum gravity, and calculate numerically the observables. After computing 
the most singular part of the observables turns out to be given by 

(tr$») = N\g c - g)l (a™ + A™ 1 +■■) (26) 

for g < g c , and 

(tr^) = N\g - g e )\ (c[ n) + <£ \ 1 +■■) (27) 

for g > g c . 

In the double scaling limit ( p6[ ) defines the puncture operator of the 2D 
quantum gravity. As long as the critical exponents in expresions (|2BJ) and 
( |2"7D are the same, we can define an analogous double scaling limit above 
the critical point given by z = N(g — g c ) J - This limit does not exist in 
the zero-dimensional matrix model, hence, it defines new effects beyond the 
original formulation of matrix model. In reference |y§ this limit is used to 
show that the stabilized model does not satisfy the nonperturbative KdV 
flow equations. 

The loops equations arise from the stabilized model as follows: We add 
to the potential W a perturbation 

W = Tr$ 2 - ^Tr$ 3 - ^Tr$ n . (28) 

3 Th 
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The potential of the stabilized theory becomes now 

V(/3 n ) = -{Tr$ 2 + # 2 Tr$ 4 - 2gTr<5> 3 + 2gTr<5> - N 
2 



+ /3 2 Tr$ 2n ~ 2 - 2/3 n Tr$" + 2g!3 n Tr<$> n+1 

+ TF E Tr$ fc Tr$"- 2 - fe } (29) 



fe=0 

now from the Hellmann-Feynmam theorem 

(f)F \ 1 n_2 
^— = -(Tr$ n )+g(Tr$ n+3 ) + ^ (T^Tr^"- 2 ^). (30) 

^P™ / /3 n =0 2 ^ fc=0 

The loop of lenght L is given by 

W(L) = Irre" (31) 
and is not difficult to prove that 

<9£ \ 



dPn+i J p n+1=0 
(32) 

then, if the energy of the ground state is zero for all values of (3 n the Hellmann- 
Feynmam theorem gives the set of Ward identities 

( dE\ 



> ^™ / /3„=0 



n = 0,...,oo (33) 



which are equivalent to the first loop equation 

V {jfi) {W{L)) = lo dJ ( W W W ( L ~ J ))' ( 34 ) 
We expect that the other loop equation can be found when we add to the 
potential W perturbation like : 

W = Tr$ 2 - ^Tr$ 3 - 2(3n ^- TJ Tr$ n * . (35) 

3 7li7l2... ■ 



erence 



The stabilized hamiltonian (|2]) is the supersymmetric hamiltonian of ref- 
171 restricted to the bosonic sector. The Ward identities of the su- 



persymetric matrix model becomes the Swchinger Dyson equation of the 
zero-dimensional matrix model. 

In the case of pure gravity the supersymmetry is broken non perturba- 
tively, then the Ward identities or the loop equation are modified by non 
perturbative corrections. 

4. Nonperturbative approach 

We come back to the question of the localization of the ground state in 
the stabilized theory, working near the critical coupling at finite values of N. 
In order to get information about the eigenvalue distribution of the ground 
state, we apply the Monte Carlo method in the computation of the path 
integral for the quantum system. In reference a Monte Carlo method was 
also applied to simulate observables of the matrix variable $. We consider 
here an alternative approach which optimizes the Monte Carlo calculation 
reducing the number of variables to the eigenvalues of $. 

The basic object to look for is the probability amplitude between two 
states and at times and T, respectively. This admits the path 



integral representation [IS 



iv 



i=l 



{ r T N 1 • 1 1 
expj-iV^ dtJ2(^K(t) + -V FP (X n (t)))\ (36) 



(*/(r)|*i(0)) = /n^(t)A(A(0))A(A(T))^ f (A(T))^(A(0)) 

n=l 

where A (A) is the Van der Monde determinant and Vpp the stabilized po- 
tential appearing in (|(|). If both and have nonvanishing projection 
over the ground state, this is the state which dominates at large T. The 
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amplitude behaves, in terms of the ground state energy E, 

(* f (T)\*m)~e- ET 



(37) 



We have performed, in practice, a discretization of the time T, such that 
tj = e j, where j — 0, • • • , N T and t Nr = eN T = T. The measure of integra- 
tion in ( j36l) becomes 

e~ S = expj-iVE E 1 (i MilhW + eV FP (A,(t J ))) 
+ X>g(A i (t )-A J -(t )) 

+ £log(A i (t JVr )-A J -(t* T ))| (38) 

This can be simulated by the Monte Carlo method, as if we were dealing 
with a two-dimensional statistical system of size N x ■ The variables are 
the eigenvalues \(tj) where % — 1, . . . , N, j = 0, . . . , A^. The kinetic term 
defines a nearest-neighbor interaction along the time axis and the determi- 
nants in ( |3T)D — the only vestige of Fermi statistics in the path integral — 
define a long-range interaction at the boundaries t = and t = T. There are 
no interactions between bulk variables along the N axis. 

The crucial point in the simulation is to reach a large enough value of 
the time T, in order to measure with confidence the ground state properties. 
We have taken steps with e = 0.1, and found that above T = 8 observables 
measured with (^) do not show significant contribution from the first excited 
states. We have imposed boundary conditions in the form of constant wave 
functions at both ends of the time interval ^i(X) = ^/(A) = const. , which 
have certainly nonvanishing overlap with the ground state. 
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We have implemented the Monte Carlo simulation in a VAX 9000 machine 
with Metropolis algorithm. We have made a random choice of the point 
U each time, updating then the N variables at that site. The number of 
iterations between measures has been 2000 MC sweeps and we have left a 
thermalization period equivalent to 5 x 10 4 MC sweeps of all the variables 
in the lattice. We have simulated systems with N = 5 and N = 10, starting 
with different initial conditions for the set of eigenvalues. In some of the 
simulations part of the eigenvalues were initially confined in the well around 
the local minimum and in the others all of them were in the central well. 
Irrespective of these different choices, after the thermalization period we 
ended up with a distribution of eigenvalues equal to zero in the well around 
the local minimum. This happened even for values of g very close to g c . 
Figure 2 shows a typical distribution for N = 5 and g — g c = 10~ 4 after 
8 x 10 6 MC sweeps of all the lattice. This supports strong evidence that at 
finite N the Fermi level is placed below the level of the local minimum of the 
potential, all the way up to g c . 

5. Conclusions 

We have shown that in the framework of the 1/N approximation there 
is only one perturbative vacuum. The Fermi level in the weak coupling 
regime is below the local minimum of the potential and, hence, instanton 
configurations which start or end at the local minimun do not exist. In the 
nonperturbative approach we have also made plausible that the ground state 
does not bear tunnelling between different wells of the potential. The expla- 
nation of this picture should be the following. The stochastic stabilization 
may be viewed as the usual stochastic quantization with asymptotic final 
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and initial states fixed and given by well denned configurations of the origi- 
nal model []19]] . Therefore, in the stabilized matrix model all configurations 
have to start and end at the main well of the stabilized potential. The Fermi 
level has to be placed below the local minimun and an instanton which starts 
or ends at the local minimun cannot exist. 

One possible way of understanding the nonperturbative effects may be 
envisaged as follows. There are two kind of static solutions of the classical 
equation in the large N limit: (a) all the fermions are restricted to the main 
well of the potential, and (b) a fermion is placed at the local minimun. These 
two solutions are degenerate at first order in 1/N. Hence we expect that the 
time dependent solution of the classical equation of motion connects the main 
well and the local minimun. The perturbative effects lift the degeneracy, and 
the solution with a fermion in the local minimun becomes an excited state. 
The nonperturbative effects have to arise from trajectories of one fermion 
which start and end at the main well. These closed trajectories are given 
by a sucession of instantons and anti-instantons |2IJ . Hence, the instantonic 
action of one fermion is 

S mst = 2 / J2(V - E F ) (39) 

J a 



which agrees with the instantonic action calculated in JT0| . 

In this paper we have considered the double scaling limit in the order: 
first take the large N limit and then g — > g c . This is the only double scaling 
limit defined in the zero-dimensional matrix model. But in the stabilized 
model it is possible to perform the double scaling limit from finite values of 
N and g — g c . We have performed a preliminary numerical calculation. In 
order to understand how an alternative double scaling limit may be achieved 
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from finite values of N, a more detailed investigation of the phases in the 
space of parameters (N,g) should be carried out. 
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Appendix A 

In this appendix we show how formulas like 

d r b 1 

" / d\ . (40) 

which have end points singular integrand, can be calculated. 
We define two fixed arbitrary points Ai and A 2 such that 

rb I r A 2 I Mi 1 



/ dX = = / dX - =+ / 



2(E y p>-V) J " ^2(E ( f>-V) JAa s]2{Ef -V) 

fb I 
+ / dA . = = (41) 



For g ^ g c , ^ has only one zero in the interval [a, b], then we choose Ai 

dV 



and A 2 such that the zero of ^ is placed between Ai and A 2 , see (fig.l, a). 
The second integral can be derived directly 

d Mi i Mi 9 f 1 ) 

^ / dX . = / dA^- ^ ; I (42) 

and the result is well defined for all values of the coupling constant. 

The other integrals are posibly singular and we perform them as follows 

fb i f b i v 

/ dX . = = / dX , . 

7ai yj2{Ef - V) J ^ ^2(EP -V) V 

= -£4^ 2(40) - y)) (43) 

where the dots are A derivatives, and integrations by parts gives 

b 1 
dX 



Ai 



1 = = Analytic terms + T dX^r ( i-} \/2(eP - V) 
v/ 2 (4 0) - V) ^ dA 



(44) 
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the result is well defined for g ^ g c and can be derived directly, but if g = g c 
then V(b) =0 and the integral is posibly singular. This is the origin of the 
critical point g c , which is defined by the condition V(b) = 0, see (fig.l,c). 

To all orders in 1/N, the energy and observables are some combinations 
of integrals which have the form 

d> d k I ft P(\) \ 

<Wd*[J. PiE-V m ) E __ El> __ a 

This integrals have the same critical point g c . The critical point g c remains 
constant to all orders in 1/N expansion. 
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Figure Captions 

Figure 1. Plot of the stabilized potential a) below the critical point, b) above the 
critical point and c) at the critical point 

Figure 2. Plot of the normalized fermionic density: continuun line. For N — 5, 
g — g c = 10 -4 , 100 time intervals . We have performed 4200 measures. 
The planar fermionic density is given by the dashed line. The vertical 
dashed line is placed at the absolute minimun of the stabilized potential. 
The local maximun and the cut of the support of p are very near of the 
local minimun. 
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